skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Calkins, Michael"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. We analyse the results of direct numerical simulations of rotating convection in spherical shell geometries with stress-free boundary conditions, which develop strong zonal flows. Both the Ekman number and the Rayleigh number are varied. We find that the asymptotic theory for rapidly rotating convection can be used to predict the Ekman number dependence of each term in the governing equations, along with the convective flow speeds and the dominant length scales. Using a balance between the Reynolds stress and the viscous stress, together with the asymptotic scaling for the convective velocity, we derive an asymptotic prediction for the scaling behaviour of the zonal flow with respect to the Ekman number, which is supported by the numerical simulations. We do not find evidence of distinct asymptotic scalings for the buoyancy and viscous forces and, in agreement with previous results from asymptotic plane layer models, we find that the ratio of the viscous force to the buoyancy force increases with Rayleigh number. Thus, viscosity remains non-negligible and we do not observe a trend towards a diffusion-free scaling behaviour within the rapidly rotating regime. 
    more » « less
  2. Abstract We present the discovery of 11 new transiting brown dwarfs (BDs) and low-mass M dwarfs from NASA’s Transiting Exoplanet Survey Satellite (TESS) mission: TOI-2844, TOI-3122, TOI-3577, TOI-3755, TOI-4462, TOI-4635, TOI-4737, TOI-4759, TOI-5240, TOI-5467, and TOI-5882. They consist of five BD companions and six very-low-mass stellar companions ranging in mass from 25MJto 128MJ. We used a combination of photometric time-series, spectroscopic, and high-resolution imaging follow-up as a part of the TESS Follow-up Observing Program (or TFOP) to characterize each system. With over 50 transiting BDs confirmed, we now have a large enough sample to directly test different formation and evolutionary scenarios. We provide a renewed perspective on the transiting “brown dwarf desert” and its role in differentiating between planetary and stellar formation mechanisms. Our analysis of the eccentricity distribution for the transiting BD sample does not support previous claims of a transition between planetary and stellar formation at ∼42MJ. We also contribute a first look into the metallicity distribution of transiting companions in the range 7–150MJ, showing that this does not support a ∼42MJtransition too. Finally, we also detect a significant lithium absorption feature in one of the BD hosts (TOI-5882). However, we determine that the host star is likely old based on rotation, kinematic, and photometric mdeasurements. We therefore claim that TOI-5882 may be a candidate for planetary engulfment. 
    more » « less
    Free, publicly-accessible full text available July 4, 2026
  3. A numerical investigation of an asymptotically reduced model for quasigeostrophic Rayleigh-Bénard convection is conducted in which the depth-averaged flows are numerically suppressed by modifying the governing equations. At the largest accessible values of the Rayleigh number Ra, the Reynolds number and Nusselt number show evidence of approaching the diffusion-free scalings of Re ∼ RaE/Pr and Nu ∼ Pr−1/2Ra3/2E2, respectively, where E is the Ekman number and Pr is the Prandtl number. For large Ra, the presence of depth-invariant flows, such as large-scale vortices, yield heat and momentum transport scalings that exceed those of the diffusion-free scaling laws. The Taylor microscale does not vary significantly with increasing Ra, whereas the integral length scale grows weakly. The computed length scales remain O(1) with respect to the linearly unstable critical wave number; we therefore conclude that these scales remain viscously controlled. We do not find a point-wise Coriolis-inertia-Archimedean (CIA) force balance in the turbulent regime; interior dynamics are instead dominated by horizontal advection (inertia), vortex stretching (Coriolis) and the vertical pressure gradient. A secondary, subdominant balance between the Archimedean buoyancy force and the viscous force occurs in the interior and the ratio of the root mean square (rms) of these two forces is found to approach unity with increasing Ra. This secondary balance is attributed to the turbulent fluid interior acting as the dominant control on the heat transport. These findings indicate that a pointwise CIA balance does not occur in the high Rayleigh number regime of quasigeostrophic convection in the plane layer geometry. Instead, simulations are characterized by what may be termed a nonlocal CIA balance in which the buoyancy force is dominant within the thermal boundary layers and is spatially separated from the interior Coriolis and inertial forces. 
    more » « less
  4. Dynamos driven by rotating convection in the plane layer geometry are investigated numerically for a range of Ekman number ( $$E$$ ), magnetic Prandtl number ( $Pm$ ) and Rayleigh number ( $Ra$ ). The primary purpose of the investigation is to compare results of the simulations with previously developed asymptotic theory that is applicable in the limit of rapid rotation. We find that all of the simulations are in the quasi-geostrophic regime in which the Coriolis and pressure gradient forces are approximately balanced at leading order, whereas all other forces, including the Lorentz force, act as perturbations. Agreement between simulation output and asymptotic scalings for the energetics, flow speeds, magnetic field amplitude and length scales is found. The transition from large-scale dynamos to small-scale dynamos is well described by the magnetic Reynolds number based on the small convective length scale, $$\widetilde {Rm}$$ , with large-scale dynamos preferred when $$\widetilde {Rm} \lesssim O(1)$$ . The magnitude of the large-scale magnetic field is observed to saturate and become approximately constant with increasing Rayleigh number. Energy spectra show that all length scales present in the flow field and the small-scale magnetic field are consistent with a scaling of $$E^{1/3}$$ , even in the turbulent regime. For a fixed value of $$E$$ , we find that the viscous dissipation length scale is approximately constant over a broad range of $Ra$ ; the ohmic dissipation length scale is approximately constant within the large-scale dynamo regime, but transitions to a $$\widetilde {Rm}^{-1/2}$$ scaling in the small-scale dynamo regime. 
    more » « less
  5. Abstract Cold Jovian planets play an important role in sculpting the dynamical environment in which inner terrestrial planets form. The core accretion model predicts that giant planets cannot form around low-mass M dwarfs, although this idea has been challenged by recent planet discoveries. Here, we investigate the occurrence rate of giant planets around low-mass (0.1–0.3M) M dwarfs. We monitor a volume-complete, inactive sample of 200 such stars located within 15 pc, collecting four high-resolution spectra of each M dwarf over six years and performing intensive follow-up monitoring of two candidate radial velocity variables. We use TRES on the 1.5 m telescope at the Fred Lawrence Whipple Observatory and CHIRON on the Cerro Tololo Inter-American Observatory 1.5 m telescope for our primary campaign, and MAROON-X on Gemini-North for high-precision follow up. We place a 95% confidence upper limit of 1.5% (68% confidence limit of 0.57%) on the occurrence ofMPsini> 1MJgiant planets out to the water snow line and provide additional constraints on the giant planet population as a function ofMPsiniand period. Beyond the snow line (100 K <Teq< 150 K), we place 95% confidence upper limits of 1.5%, 1.7%, and 4.4% (68% confidence limits of 0.58%, 0.66%, and 1.7%) for 3MJ<MPsini< 10MJ, 0.8MJ<MPsini< 3MJ, and 0.3MJ<MPsini< 0.8MJgiant planets, respectively; i.e., Jupiter analogs are rare around low-mass M dwarfs. In contrast, surveys of Sun-like stars have found that their giant planets are most common at these Jupiter-like instellations. 
    more » « less
  6. We report the discovery of TOI-4641b, a warm Jupiter transiting a rapidly rotating F-type star with a stellar effective temperature of 6560 K. The planet has a radius of 0.73 RJup, a mass smaller than 3.87 MJup(3σ), and a period of 22.09 d. It is orbiting a bright star (V=7.5 mag) on a circular orbit with a radius and mass of 1.73 R⊙ and 1.41 M⊙. Follow-up ground-based photometry was obtained using the Tierras Observatory. Two transits were also observed with the Tillinghast Reflector Echelle Spectrograph, revealing the star to have a low projected spin-orbit angle (λ=$$1.41^{+0.76}_{-0.76}$$°). Such obliquity measurements for stars with warm Jupiters are relatively few, and may shed light on the formation of warm Jupiters. Among the known planets orbiting hot and rapidly rotating stars, TOI-4641b is one of the longest period planets to be thoroughly characterized. Unlike hot Jupiters around hot stars which are more often misaligned, the warm Jupiter TOI-4641b is found in a well-aligned orbit. Future exploration of this parameter space can add one more dimension to the star–planet orbital obliquity distribution that has been well sampled for hot Jupiters. 
    more » « less
  7. ABSTRACT Numerical simulations are used to investigate large-scale (mean) magnetic field generation in rotating spherical dynamos. Beyond a certain threshold, we find that the magnitude of the mean magnetic field becomes nearly independent of the system rotation rate and buoyancy forcing. The analysis suggests that this saturation arises from the Malkus-Proctor mechanism in which a Coriolis-Lorentz force balance is achieved in the zonal component of the mean momentum equation. When based on the large-scale magnetic field, the Elsasser number is near unity in the saturated regime. The results show that the large and small magnetic field saturate via distinct mechanisms in rapidly rotating dynamos, and that only the axisymmetric component of the magnetic field appears to follow an Elsasser number scaling. 
    more » « less